GLOBAL ATTRACTOR OF 
A COUPLED TWO-CELL BRUSSELATOR MODEL 
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Dedicated to George R. Sell on the Occasion of His 70th Birthday 

Abstract. In this work the existence of a global attractor for the solution semiflow 
of the coupled two-cell Brusselator model equations is proved. A grouping estimation 
method and a new decomposition approach are introduced to deal with the challenge 
in proving the absorbing property and the asymptotic compactness of this type of four- 
variable reaction-diffusion systems with cubic autocatalytic nonlinearity and with linear 
coupling. It is also proved that the Hausdorff dimension and the fractal dimension of the 
global attractor are finite. 



1. Introduction 

Consider a coupled two-cell model of reaction-diffusion systems with Brusselator kinetics 
[12,18,20,31], 

— = diAu + a- (b+l)u + u^v + Diiw - u), (1.1) 
ot 

Qd 

— = d2Av + bu-u^v + D2(z-v), (1.2) 
ot 

— = diAw + a- {b+l)w + w^z + Di{u - w), (1.3) 
d z 

— = d2Az + bw - w'^z + D2(v - z), (1.4) 
ot 

for t > 0, on a bounded domain C 5R",n < 3, that has a locally Lipschitz continuous 
boundary, with the homogeneous Dirichlet boundary condition 

u{t,x) = v{t,x) = w{t,x) = z{t,x) = 0, t> 0, X € dft, (1.5) 

and an initial condition 

n(0, x) = no(x), ^(O, x) = vo{x), w{0, x) = wo{x), z{0, x) = zo{x), x £ (1-6) 

where di,d2,a,b, Di, and D2 are positive constants. We do not assume that the initial 
data uo,vo,wo, zq, nor the solutions u{t,x),v{t,x),w{t,x), z{t,x) to be nonnegative. In 



2000 Mathematics Subject Classification. 37L30, 35B40, 35B41, 35K55, 35K57, 35Q80, 80A32, 92B05. 
Key words and phrases. Brusselator, two-cell model, global dynamics, global attractor, absorbing set, 
asymptotic compactness. 

1 



2 



YUNCHENG YOU 



this work, we shall study the asymptotic dynamics of the solution semiflow generated by 
this problem. 

The Brusselator is originally a system of two ordinary differential equations as a model 
for cubic autocatalytic chemical or biochemical reactions, cf. [2,25,35]. The name is after 
the hometown of scientists who proposed it. Brusselator kinetics describes the following 
scheme of chemical reactions 

A — ^ U, 
B + U — > V + D, 
2U + V — > 3U, 
U — ^ E, 

where A, B, D, E, U, and V are chemical reactants or products. Let u{t,x) and v{t,x) be 
the concentrations of U and V, and assume that the concentrations of the input compounds 
A and B are held constant during the reaction process, denoted by a and b respectively. 
Then one obtains a system of two nonlinear reaction-diffusion equations called (diffusive) 
Brusselator equations, 

^ = diAu + u^v -{b+l)u + a, (1.7) 

^ = d2Av - u^v + bu, (1.8) 
ot 

There are several known examples of autocatalysis which can be modeled by the Brusse- 
lator equations, such as ferrocyanide-iodate-sulphite reaction, chlorite-iodide-malonic acid 
reaction, arsenite-iodate reaction, some enzyme catalytic reactions, and fungal mycelia 
growth, cf. [1,2,5,9]. 

Numerous studies by numerical simulations or by mathematical analysis, especially after 
the publications [22,23] in 1993, have shown that the autocatalytic reaction-diffusion sys- 
tems such as the Brusselator equations and the Gray-Scott equations [13, 14] exhibit rich 
spatial patterns (including but not restricted to Turing patterns) and complex bifurca- 
tions [1,4,5,7,8,11,16,24,26,38] as weh as interesting dynamics [6,9,10,17,21,27-29,39,40] 
on ID or 2D domains. 

For Brusselator equations and the other cubic autocatalytic model equations of space 
dimension n < 3, however, we have not seen substantial research results in the front 
of global dynamics until recently this author proved the existence of a global attractor 
for Bruuselator equations [41], Gray-Scott equations [42], Selkov equations [43], and the 
reversible Schnackenberg equations [44]. 

In this paper, we shall show the existence of a global attractor in the product phase 
space for the solution semiflow of the coupled two-cell Brusselator equations (ll.iP " (ll.4p 
with homogeneous Dirichlet boundary conditions (II. 5p . 

This study of global dynamics of the two-cell model of four coupled components is a 
substantial advance from the one-cell model of two-component reaction-diffusion systems 
toward the biological network dynamics [12,19]. Multi-cell models generically mean the 
coupled ODEs or PDEs with large number of unknowns (components), which appear widely 
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in the literature of systems biology as well as cell biology. Here understandably "cell" is a 
generic term that may not be narrowly or directly interpreted as a biological cell. Coupled 
cells with diffusive reaction and mutual mass exchange are often adopted as model systems 
for description of processes in living cells and tissues, or in distributed chemical reactions 
and transport for compartmental reactors [31,37]. The mathematical analysis combined 
with semi-analytical simulations seems to become a common approach to understanding 
the complicated molecular interactions and signaling pathways in many cases. 

In this regard, unfortunately, the problem with high dimensionality can arise and puzzle 
the research when the number of molecular species in the system turns out to be very 
large, which makes the behavior simulation extremely difficult or computationally too 
expensive. Thus theoretical results on multi-cell model dynamics can give insights to 
deeper exploration of various signal transductions and tempro-spatial pattern formations. 

For most reaction-diffusion systems consisting of two or more equations arising from the 
scenarios of autocatalytic chemical reactions or biochemical activator-inhibitor reactions, 
such as the Brusselator equations and the coupled two-cell Brusselator equations here, the 
asymptotically dissipative sign condition in vector version 

lim F(s) • s < 

is inherently not satisfied by the opposite-signed and coupled nonlinear terms, see (jl.lip 
later. Besides a serious challenge in dealing with this coupled two-cell model is that, 
due to the coupling of the two groups of variables n, v and w, z, one can no longer make 
a dissipative a priori estimate on the w-component by using the iJ-equation separately 
and then use the sum y{t, x) = u{t, x) + v{t, x) separately to estimate the u-component 
in proving absorbing property and in proving asymptotical compactness of the solution 
semiflow as we did in [41-43]. The novel mathematical feature in this paper is to overcome 
this coupling obstacle and make the a priori estimates by a method of grouping estimation 
combined with a new decomposition approach. 

We start with the formulation of an evolutionary equation associated with the two-cell 
Brusselator equations. Define the product Hilbert spaces as follows, 

H = L^iSl) X L^iyt) X L^{^) X L^{n), 

E = H^{n) X H^{n) X H^{n) x H^{n), 
n = {H^{n) n H\n)) x {H^{n) n H\n)) x {H^{n) n H^{n)) x {H^{n) n H^{n)). 

The norm and inner-product of H or the component space L'^(Q) will be denoted by || • || 
and (•,•), respectively. The norm of L'f{i}) will be denoted by || • \\lp if p ^ 2. By the 
Poincare inequality and the homogeneous Dirichlet boundary condition (jl.Sp . there is a 
constant 7 > such that 

\\y^f>l\\^f, foTipe H^{n) or E, (1.9) 

and we shall take || V(/?|| for the equivalent norm of the space E and of the component space 
Hq{Q). We use | • | to denote an absolute value or a vector norm in a Euclidean space. 
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It is easy to check that, by the Lumer-PhiUips theorem and the analytic semigroup 
generation theorem [34], the hnear operator 
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V 












diA 
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:DiA)i= U)^H 



(1.10) 



is the generator of an analytic Co-semigroup on the Hilbert space H, which will be denoted 
by e^\t > 0. By the fact that H^{n) L^{Q) is a continuous embedding for n < 3 and 
using the generalized Holder inequality. 



\w 
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^ II^IIls Ihllis, \\w z\\ < 
one can verify that the nonlinear mapping 

/a - {b + l)u + u^v + Di{w 
bu — v?v + D2{z — v) 



zWia^ for n, -w, z G L^(r2), 



F{g) 



u)\ 



w) 



\ 



: E 



(1.11) 



[b+ l)w + w'^z + Di{u 
bw — nP'z + D2{v — z) J 

where g = {u,v,w,z), is well defined on E and locally Lipschitz continuous. Then the 
initial-boundary value problem ()l.ip - (11.6p is formulated into the following initial value 
problem, 

^ = Ag + F{g), t > 0, (1.12) 

5(0) = go = col {uo,Vo,Wo, Zo). 
where g{t) = col (u(t, ■),v{t, •),w{t, •), z{t, •)), simply written as (u(t, •),v{t, ■),w{t, •), z{t, •)). 
Accordingly we shall write go = col {uq, vo,wq, zq). 

By conducting a priori estimates on the Galerkin approximate solutions of the initial 
value problem (jl.l2p and the weak convergence, we can prove the local existence and 
uniqueness of the weak solution g{t) of (jl.l2p in the sense of J. M. Ball specified in [3], 
which then is shown to be a local mild solution [3] and further turns out to be a local strong 
solution [34, Theorem 46.2]. Moreover, by taking the i?-inner-product of (11.121) with this 
strong solution g(t) itself and conducting a priori estimates, one can prove the continuous 
dependence of the solutions on the initial data and the following property, which is satisfied 
by the strong solution g, 

g G C{[0, Tmax)\H) n Ci((0, Tmax);H) D L\0, Tmax;E), (1.13) 

where Imax = [0,Tjnax) is the maximal interval of existence. 

We refer to [15, 34, 36] and many references therein for the concepts and basic facts in 
the theory of infinite dimensional dynamical systems, including few given below for clarity. 

Definition 1. Let {S'(t)}t>o be a semifiow on a Banach space X. A bounded subset Bq 
of X is called an absorbing set in X if, for any bounded subset C X, there is some finite 
time to > depending on B such that S{t)B C Bq for all t > to. 
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Definition 2. Let {S{t)}t>o be a semiflow on a Banach space X. A subset £/ of X is called 
a global attractor for this semiflow, if the following conditions are satisfied: 

(i) £/ is a nonempty, compact, and invariant set in the sense that 

S{t)£/ = £/ for any t > 0. 

(ii) attracts any bounded set i? of X in terms of the Hausdorff distance, i.e. 

dist{S{t)B, .s/) = sup inf \\x — y\\x — > 0, as i — > oo. 

Definition 3. A semiflow {S{t)}t>o on a Banach space X is called asymptotically compact 
if for any bounded sequences {xn} in X and {tn} C (0,oo) with t„ oo, there exist 
subsequences {xn,.} of {un} and {tn^} of {tn}, such that linif^^rx S{tn,.)xnk exists in X. 

Here is the main result of this paper. We emphasize that this result is established 
unconditionally, neither assuming initial data or solutions are nonnegative, nor imposing 
any restriction on any positive parameters involved in the equations (jl.ip - (jl.4p . 

Theorem 1 (Main Theorem). For any positive parameters di,d2,a,b,Di,and D2, there 
exists a global attractor in the phase space H for the solution semiflow {S{t)}t>o gen- 
erated by the coupled two-cell Brusselator evolutionary equation (11.12p . 

For investigation of the asymptotic compactness for the two-dell Brusselator semiflow, 
we shall take the approach of showing the ^-contracting property. Recall the definition of 
the Kuratowski measure of noncompactness for bounded sets in a Banach space X, 

k{B) '= inf {6 : B has a finite cover by open sets in X of diameters < (5} . 

If B is an unbounded set, then we define k,{B) = 00. The basic properties of the Kuratowski 
measure are listed here, cf. [34, Lemma 22.2]. 

Let X be a Banach space and k be the Kuratowski measure of noncompactness of 
bounded sets in X. Then k has the following properties: 

(i) k{B) = if and only if B is precompact in X, i.e. ClxB is a compact set in X. 

(ii) k{Bi) < k{B2) whenever Bi C B2. 

(ili) k{Bi + B2) < k{Bi) + k{B2), for any linear sum Bi + B2- 

The following lemma states concisely the basic result on the existence of a global at- 
tractor for a semiflow and provides the connection of the K-contracting concept to the 
asymptotical compactness, cf. [34, Chapter 2]. 

Lemma 1. Let {S{t)}t>o be a semiflow on a Banach space X. // the following conditions 
are satisfied: 

(i) {S'(t)}j>o has a bounded absorbing set in X, and 

(ii) {S{t)}t>o is K- contracting, i.e. limt_,oo K,{S{t)B) = for any bounded set B G X, 
then {S{t)}t>o is asymptotically compact and there exists a global attractor £/ inX for this 
semiflow. The global attractor is given by 

^ = iv{Bo) = f] Clx(j{S{t)Bo). 

r>0 t>T 
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In Section 2 we shall prove the global existence of the strong solutions of the two- 
cell Brusselator evolutionary equation (jl.l2p and the absorbing property of this coupled 
Brusselator semiflow. In Section 3 a new decomposition technique is presented to deal with 
the asymptotic compactness issue of this problem. This approach is taken to show the k- 
contracting property for the {v, z) components in Section 4 and for the (n, w) components 
in Section 5, respectively. In Section 6 we assemble these results to prove the existence of a 
global attractor in the phase space H for the coupled Busselator semiflow and to show that 
the global attractor has a finite Hausdorff dimension and a finite fractal dimensions. As a 
remark, with some adjustment in proof, these results are also valid for the homogeneous 
Neumann boundary condition. Furthermore, corresponding results are valid for the coupled 
two-cell Gray-Scott equations, Selkov equations, and Schnackenberg equations. 

2. Absorbing Property 

In this paper, we shall write u{t,x),v{t,x),w{t,x), and z{t,x) simply as u(t),v{t),w{t), 
and z{t), or even as u,v,w, and z, and similarly for other functions of {t,x). 

Lemma 2. For any initial data go = {uo,vo,wo, zq) G H, there exists a unique, global, 
strong solution g{t) = {u{t),v{t),w{t),z{t)), t € [0,oo), of the coupled Brusselator evolu- 
tionary equation (I1.12p . 

Proof. Taking the inner products ( <\1.2^ .v(t)) and ( (|1.4p . z(t)) and summing them up, we 
get 

U^Avf + i\\4']+d2(\\Vv\\' + \\Vvf) 



2 \dt" " dt' 



J (— n V + buv — w z + bwz — D2[v — 2vz + z ]) dx ^2 1) 



n 



dx + -b'^\n\ < -b'^ini. 
2 ' ' - 2 ' ' 



+ \\zf) +2jd2{\\vf + \\zf) <b^\n\ 



It follows that 

d_ 

di 

which yields 

Mt)f + \\z{t)f < e-^-rd^t (||^^||2 ^ ||^^||2) ^ f^^^ ^ [0,Tma.). (2.2) 

Let y{t, x) = u{t, x) + v{t, x) + w{t, x) + z{t, x). In order to treat the u-component and 
the i(;-component, first we add up (jl.ip . (|1.2p . (jl.3p and (jl.4p altogether to get the following 
equation satisfied by y{t) = y{t,x), 

^ = diAy -y + [{d2 - di)A{v + z) + {v + z) + 2a] . (2.3) 
Taking the inner-product ( (|2.3p . -»(t)) we obtain 
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~\\y\? + dillVyf + \\yf = [{d2 - d,)Aiv + z) + {v + z) + 2a]ydx 

< \di - d2\\\V{v + z)\\\\Vy\\ + \\v + z\\\\y\\ + 2a\n\'/^\\y\\ 

< y livyf + ^-^^^IMv + z)f + l\\yf + \\v + zf + 4a'\n\, 

so that 

+ djV.f + \\yf < ^^^l|V(. + + 4 + ) + 8a'\n\. 

Then we get 

+ dillVyf + < '^^~^^'V (^ + + C^i(^o, ^o,t), (2.4) 

where 

Ci(t;o, ^o,t) = 4e-2^'^^* {\\vof + \\zof) + + 8a') M- (2.5) 

Integrate the inequahty ()2.4p to see that the strong solution y(t) of (j2.3p satisfies the 
following estimate, 



Mt)f < \\uo + VO + WO + zof + ^^^^ r l|V(Ks) + ^(5))f ds 







2 / 46^ 



(2.6) 



From (12. ip we have 

d2 f \\V{v{s) + z{s))f ds < 2d2 f {\\Vv{s)f + \\Vz{s)f) ds < {\\vof + \\zof) +b''\n\t. 
Jo Jo 



Substitute this into (12.61) to obtain 

4_ _|- I — 

di d2 7(^2 



'<\\uo + vo + wo + zof+( + ^ ) (Ibof + "~ "'^ 



A "i"2 702 / J 
Let p(t) = u{t) + tt;(t). Then by (|2.2p and (j2.7p we have shown that 

= \\u{t) + w{t)f = \\y{t) - {v{t) + z{t))f 



(2.7) 



(2 8) 

<2{\\uo + vo + wo + zof + \\vo\\^ + \\zof) +C2{go)t, foite[0,Tmax), 
where C2{go) is a constant depending on the initial data qq. 
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On the other hand, let tp{t,x) = u(t,x) + v{t,x) — w{t,x) — z{t,x), which satisfies the 
equation 

^ = diA^ - (1 + 2Di)iP + [{d2 - di)A{v -z) + {l + 2{Di - D2)){v - z)] . (2.9) 
Taking the inner-product ( (j2.9p . -^(t)) we obtain 

~Uf + rfillV^f + ll^f < ~Uf + dillV^f + (1 + 2Dr)Uf 
< {di - d2)\\V{v - z)\\\\V^\\ + |1 + 2(Z)i - D2)\\\v - z\\m 



so that 



where 



+di||VV^f + llV^f < ^l^-^||V(t; - z)f + C3(t;o,^o,i), (2.10) 



C,{v^, zo,t) = 2|1 + 2iD, - D2)\^ (e-^^-^^* {Hf + \\zof) + 
Integration of (|2.10p yields 

2 < II^Q + - - Zof + ^1^^ /* ||V(t-(s) - Z{s))f ds 



(2.11) 



1 Jo 



|i + 2(Di - Z)2)r ( ^(Ibof + Ikof ) + ^ i ) ' t e [o,r„,,). 



(2.12) 



Note that 
ft 



d2 [ \\Viv{s)-z{s))fds<2d2 [ {\\Vv{s)f + \\Vz{s)f) ds<{\\vof + \\zof)+b''\n\t. 
Jo Jo 

From ([21^ it follows that 

< \\uo + VO-WO- zof + 'VJ^''' dl^oll' + ll^oll' + ^'1^1 

+ |1 + 2iDi - D2)\' f ^(ll^of + \\zof) + ^-^t), te [0,T^a.). 
\jd2 7^2 / 

Let q{t) = u{t) - w{t). Then by (g^l) and (I2T3I1 we find that 

\\q{t)f = Mt) - w{t)f = \\m - iv{t) - z{t))f 

<2{\\uo + vo-wo- zof + \\vof + \\zof) +Ci{go)t, for t G [0, T^^x), 

where Ci{go) is a constant depending on the initial data qq. 

Finally combining (12.81) and (I2.14P we can conclude that for each initial data go € H, 
both u{t) = {l/2){p{t) + q{t)) and w{t) = (l/2)(p(t) — q{t)) components are bounded if 
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Tmax of the maximal interval of existence of the solution is finite. Together with (j2.2p . 
this shows that, for each qq £ H, the strong solution g{t) = {u{t),v(t),'w{t), z{t)) of the 
equation ()1.12p will never blow up in H at any finite time and it exists globally . □ 

Due to Lemma O the family of all the global strong solutions {g{t;go),t > 0,go £ H} 
defines a semiflow on H, 

S{t) : go^ g{t;wo), goeH,t>0, 

which is called the coupled Brusselator semiflow generated by the coupled Brusselator 
evolutionary equations. 

Lemma 3. There exists a constant Kq > 0, such that the set 

Bo = {\\g\\ € H : \\gf < Ko} (2.15) 
is a bounded absorbing set Bq in H for the coupled Brusselator semiflow {S{t)}t>o- 
Proof. For this coupled Brusselator semiflow, from (|2.2p we obtain 

limsup(||T;(t)f + \\z{t)f) <Ro = (2.16) 



Moreover, for any t > 0, (j2.ip also implies that 

'^\\\Vv{s)f + \\Vz{s)f) ds < U\Ht)f + \\z{t)f + h^M 

«2 



1 f 2'ydofn, ,,2 n n2N b'^l^W b'^l^l 

< :r e^^dl^of + ll^of) + 



d2 V 27(i2 



(2.17) 



which is for later use. 

From (12.41) we can deduce that 



4 (e*l|y(i)f) < ^^-^ e*||V(^;(t) + z{t))f + e'Ci{vo, zo,t). (2.18) 



dt ' d 

Integrate (12.181) to obtain 

\\y{t)f <e~*||'Uo + wo + u'o + ^o 



+ ^1-^ /* e-(*--) II V(t;(r) + z(T))f dr + C,{vo, 
"1 Jo 



(2.19) 



where 



C5{vo,zo,t) = e-' l\e^'-^'"''^-dT{\\vof + \\zof) + (^^+80^^ \n\ 
<4a(0(|kof + lkof)+ (^+8a2) |0|, 
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a 



f ^4^^"'"'^*' if 1 - 27^2 > 0; 
te-* < 2e-ie'*/2, if 1 - 27^3 = 0; 

if 1 - 27^2 < 0. 



(2.20) 



' re-*, 



11-27^2! 

On the other hand, multiplying (j2.ip by e* and then integrating each term of the resulting 
inequality, we get 



. d 



1. 



v{rW + MT)f) dr + d, e^{\\Vv{T)f + \\z{T)f)dT<-bMe\ 



so that, by integration by parts and using (j2.2p . we obtain 



1 



d2 / e{\\Vv{T)r + \\yv{T)r)dT<-hMe'-- I e-[\^v{T)f + \\Vz{T)f)dT 



d 



e\\\v{t)f + \\z{t)f 



+ 



<h-'\n\e' + {\\vof + 



+ 



,(l-27d2)T 



+ 



e^{\\v{T)f + \\z{T)f)dT 



) dr + ^r-^e* 



27^2 



< 1 + 



27^2 



62|0|e* + (l + a(t)e*)(||^;of + 



for t > 0. 



(2.21) 



Substituting (p:2T]) into ([2T9]) . we obtain that, for t > 0, 
'2 < e-*||'Uo + t'o + 1L'o + 2:o|P + C'5('yo,^0,i) 



+ 



2|di-(i2p ^-t 

di d2 



1 + 



1 



27^2 



b^\n\e' + {l + e'ait))i\\vof + 



< e *||uo + 7;o + -wo + zof + 4a(t)(||t;o||^ 



V(f + 8a^l|fi| 



(2.22) 



I 2\d,-d2\^ ^_t 

di d2 



1 + 



27(^2 



62|17|e* + (l + e*a(t))(||t;of + 



Note that (121201) shows a{t) ^ 0, as t ^ 0. From (121221) we find that 



46^ 

limsup ||y(t)f < i?i = 1 + ( — + 8a^ ) + 



t— >oo 



7^2 



2|(ii -dal 
di d2 



1 + 



27^2 



(2.23) 



The combination of (|2.16|) and (|2.23|) gives us 



limsup \ \u{t) + w 

t—*oo 



limsup \\y{t) - {v{t) + z{t))f < 4i?o + 2Ri. (2.24) 

t—*oo 
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Similarly, from the inequality (j2.10p satisfied by tp(t) = u{t) + v{t) — w{t) — z{t), we get 



d 



< 



\di - d 



'2 1 t 



dt ^ ' d\ 

Integrate ()2.25p to obtain 

< e'^lllio + Wo - w^o - 2:o|P 



e*||V(t;(t)-z(t))f + e*C3(^o,^o,t). 



+ 



^(*-) ||V(r;(r) - z(r))f dr + C76(t;o, zo, t), 



(2.25) 



(2.26) 



where 



C6(T;o,zo,t) = 2|1 + 2(I)i - D^t (e-*^*e(i-2^'^^)^dT {\\v^f + 

< 2|1 + 2(Z?i - D^)\^ (a(t)(|bof + Uf) + 

Using (j2.2ip to treat the integral term in (j2.26p . we obtain that 

< e"*||no -^vq-wq -^oIP + Q (-yo , , *) 



+ 



id; 



d\ d2 



e^{\\Vv{T)f + \\Vz{T)f)dT 



< e-'\\uo + VO-WO- + 2|1 + 2{Di - Z^s)!' aWdkoIr + 



7"2 



+ 



2|di-(i2p ^_t 



1 + 



1 



27(i2 



b'\Q\e' + {l + e'ait)){\\vof + \\zof) 



Therefore, since a{t) — > 0, as t 0, from (j2.27p we get 



limsup ||V(i)f < i?2 = 1 + 2b^\n\ 

t— >oo 



|l + 2(Z)i-Z)2)P , \di-d2 



+ 



1 + 



27(i2 



7^2 rfl C?2 

The combination of (I2.16P and (I2.28P gives us 

limsup \\u{t) - w{t)f = limsup \\4^{t) - {v{t) - z{t))f < ARq + 2R2. 

t— >oo t— ►oo 

Finally, putting together ()2.24p and (|2.29p . we assert that 

limsup(||u(t)f + \\w{t)f) < 8R0 + 2{Ri + R2). 

t— >oo 

Then assembling (|2.16p and (|2.30p . we end up with 



for t > 0. 
(2.27) 

(2.28) 

(2.29) 

(2.30) 



limsup ||(7(t)|P = limsup(||u 



+ vltW + w 



+ z 



<9Ro + 2iRi+R2). 



t— >oo t— >oo 

Thus this lemma is proved with Kq = 9Rq + 2{Ri + R2) in the absorbing ball Bq in (j2.15p . 
And Kq is a uniform positive constant independent of initial data. □ 
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3. A New Decomposition for Asymptotic Compactness 

The lack of inherent dissipativity and the cross-cell coupling make the attempt of showing 
the asymptotic compactness of the coupled Brusselator semiflow even more challenging. A 
generic and good idea in dealing with this issue is through a decomposition approach. The 
existed decomposition methods in different scenarios have been commented in Section 3 
of [41] and [42]. 

In [41,42] the author proved the following lemma, which provides an effective decom- 
position approach and has been used in proving the ^-contracting property linked to the 
existence of a global attractor for several cubic autocatalytic reaction-diffusion systems 
in [41-44]. 

Lemma 4. Let {^{t)}t>o be the solution semiflow generated by the Brusselator equations 
(jl.7p - ()1.8p on = X Then there exists a global attractor £/ in "K for this 

semiflow if and only if the following two conditions are satisfied : 

(i) There exists a bounded absorbing set "Bq in 'K for this semiflow. 

(ii) For any e > 0, there are positive constants M = M{e) and T = T{e) such that 

I \S{t)wQ\^ dx <Ce, for any t>T, Woe 'Bo, (3.1) 

J Q.{\v{t)\>M) 
where C > is a uniform constant, and 

K ((5(t)'Bo)n(it,(t)|<M)) — >0, ast^oo, (3.2) 

where 

def 

(5(t)So)n(|^(t)|<Af) = {{S{t)wo){-)CM{-;t,wo) : for Wo G So} , 
in which (Mix;t,wo) is the characteristic function of the subset ^l{\v{t)\ < M), and v(t) = 
v{t,x,wo) is the v-component of the solution of the Brusselator equations (11.7p - (ll.8p . 

However, this lemma does not work for the investigation of the ^-contracting property 
with regard to the coupled Brusselator system (jl.ip - (jl.4p . 

Motivated by LemmaHJ we now prove a new decomposition technique in the next lemma, 
which relaxes the decomposing criterion depending on the truncation of one component 
function (say, the v-component) to allowing the criteria depending on the truncation of two 
component functions (say, the T;-component and the z-component), one for each cell, and 
the decomposition of the two subgroups of components can be different. But unlike Lemma 
m the next lemma is only a sufficient condition for the existence of a global attractor. 

Lemma 5. For the solution semiflow {S{t)}t>o generated by the coupled Brusselator evo- 
lutionary equation (jl.l2p on H , there exists a global attractor in H , if the following two 
conditions are satisfied: 

(i) There exists a bounded absorbing set Bo in H for this semiflow. 

(ii) For any e > 0, there are positive constants M = M{e) and T = T{e) such that 

I \v{t)\^dx+ / \z{t)\^ dx < Lie, for any t > T, go e Bo, (3.3) 

J Q.{\v{t)\>M) Jn{\z{t)\>A'I) 
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and 

I \u{t)\'^dx+ \w{t)\'^ dx < L2e, for any t > T, go e Bq, (3.4) 

Jn(\v{t)\>M) Jn{\z{t)\>M) 

where Li > and L2 > are two uniform positive constants. 
(iii) For any given M > 0, 

K {Pu,v [('S'(t)-Bo)n(|^(t)|<M)]) — >0, as t^ 00, (3.5) 

and 

K {Pw,z [(.S{t)Bo)nQ^(_t)\<M)]) — ^ 0, as t^ 00, (3.6) 

where P^^z (^nd Pu,w O'T^ respectively the orthogonal projections from H onto the component 
spaces L^(i7)^ x L'^{il,)z and L'^{^l)u x L^(r2)u,, 

(5(t)So)n(|i,(t)|<M) {{S{t)go){-)eM{-;t,go) : for 5-0 G -Bo}, (3.7) 
(S(t)i?o)n(|.(i)|<M) ='{(^(i)5o)(-)CM(-;t,5o) :for5oG-Bo}, (3.8) 

in which 9M{x;t, go) and (,M{x;t, go) are respectively the characteristic functions of the 
subsets n{\v{t)\ < M) and n{\z{t)\ < M), and v{t) = v{t,x,go) is the v-component, 
z{t) = z{t,x,gQ) is the z-component of the solutions of the coupled Brusselator evolutionary 
equations ()1.12p . 

Proof. In light of Lemma [H it suffices to show that this solution semiflow {S'(t)}t>o is 
K-contracting on the space H. Since the absorbing set Bq in (j2.15p attracts every bounded 
set B C H, we need only to show 

lim K{S{t)Bo) = 0. (3.9) 

By the linear sum property of the Kuratowski measure listed in Section 1, we have 

fiiSit)Bo) < K{Siit)Bo) + K{S2it)Bo), t > 0, (3.10) 

where 

and it holds that K{Si{t)Bo) = K{Pu,y[S{t)Bo]) and K{S2{t)Bo) = K{Pw^z[S{t)Bo]). 
Note that for any given constant M > 0, we have 

Si{t)Bo c {Si{t)Bo)eM + {Si{t)Bo){l - 9m), (3.11) 

where 

{Si{t)Bo)eM = {{S{t)go){-)eM{-,t,go) : go G Bo}, 
{Si{t)Bo)il - Om) = {{S{t)go){-){l - eM{-,t,go)) : 9o G Bo}. 
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By (|3.3p and (j3.4p . for an arbitrarily given e > 0, there exist constants M > and T > 
such that 



\{Si{t)go){x){l - eM{x;t,go))\Ux = [ \Si{t)g, 

Jn{\v{t)\>M) 

= + < Le, t>T, 

Jn( 



'n(\v{t)\>M) 

where L = Li + L2, which imphes that 

K((5i(t)So)(l-0A/) <2\/Li, t>T. (3.12) 

On the other hand, by (|3.5p . for the same e and M, there exists a sufficiently large > 0, 
such that 

K{{Si{t)Bo)eM) = k{Pu,. [(5(t)5o)n(|.(t)|<M)]) < t> T\ (3.13) 

Then by (j3.1ip and the monotone property of the K-measure, ()3.12p and (j3.13p show that 

K{Si{t)Ba)) <K{{Si{t)Bo)9M) + K((5i(t)5o)(l - Om)) <e + 2Vre, fort > max{r,T^}. 

Similarly from (13. 3p . ()3.4p . and (13. 6p we can get a sufficiently large > 0, such that 

K{S2{t)B^))<ti{{S2{t)Bo)iM) + ^{{Si{t)Bo){l-iM))<e + 2^/Te, for t > max{r, T^}. 

Finally we substitute the last two inequalities into (j3.10p to conclude that (j3.9p is valid. □ 

In the next two sections we shall check the conditions specified in the items (ii) and 
(iii) of Lemma [5] toward the proof of the existence of a global attractor for the coupled 
Brusselator semiflow. 

4. K-Contracting Property for the (u, z) Components 

In this section, we shall check that the conditions specified in the items (ii) and (iii) of 
Lemma [5] for the {v^z) components of the coupled Brusselator equations. 
In this section and next section, we shall use the notation 

0^ = 0((/)(t) >M) = {xen: <l){t, x) > M} 

0|^|,M = ^Mt)\ <M) = {xen: |0(t,x)| <M} 

where (j){t, x) is any measurable function on Q for each given t >0. If a function p{x),x £ $7, 
is in H or L^(r2), then we shall use the following norm notation 

Ilp|lo0 = / \p{x)\'^dx and \\p\\l = \p{x)\^ dx. 

A/ J n((i,{t)>M) Jn{\(f>{t)\<M) 

We can use m{S) or \S\ to denote the Lebesgue measure of a measurable subset S in Q. 
For any measurable defined on Q, let 

M, ifv9(x)>M, ^ , , fvix) + M, if^{x)<-M, 

{ip — M)+ = { and (w+M)_ = < 

^" i{^{x)<M- lo, if(^(x)>-M. 
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As a preliminary remark, since Bq in Lemma [3] is a bounded absorbing set in H for the 
coupled Brusselator semiflow {-S'(t)}f>o, there exists a constants Tq > 0, such that 

\\S{t)gof < Kq, for any t > Tq, go = {uo,vo,wo, zq) G Bq, (4.2) 
where Kq is the constant given in ()2.15p . Let this Tq be fixed. 

Lemma 6. For any e > 0, there exist positive constants Mi = Mi(e) and Ti = Ti(e), 
such that the v-component v{t) = v{t,x,gQ) and the z-component z{t) = z{t,x,gQ) of the 
solutions of the coupled Brusselator equations (11.11) - (ll.4p satisfy the following estimate, 



n{\v{t)\>Mi) 



\v{t)\'^dx + 



\z{t)\^ dx < for t > Ti, 50 G ^0, (4.3) 

^(\z{t)\>lvh) 7"2 



where Li (46^)/ (7^2) is a uniform constant. 

Proof. By (OI) . for any gQ E Bq and any t > Tq, we have ||u(t)p + \\z{t)\\'^ < Kq. Hence 
we have 

[m{Q{\v{t)\ > M)) + m{Q{\z{t)\ > M))] 



< 



n{\v(t)\>M) 



\v{t)\'^dx + 



n{\z(t)\>M) 



\v{t)\'^dx < Kq, 



so that there exists an M = M(e) > such that for any t > TQ,gQ £ Bq, 

mm<t)\ > M)) < ^ < and m(17(|z(t)| > M)) < ^ < |. (4.4) 

Taking the inner-product ( (jl.2p . {v{t) — M)_|_), where M is given in ()4.4p . we obtain 

l±\\(v-MUf + d2 [ \Viv-MU\Ux 
2 dt Jqv^ 

= - u'^v{v - M)+dx + / bu{v - M)+dx + D2 {z - v){v - M)+ dx 



< 



u{v - M)+ 



+ u^M{v - M)+ 



"M 
k2 



dx + —m{n{v{t) > M)) 



+ D2 I {z-M){v- M)+ dx + D2M I 
-D2 I {v- M)l dx - D2M I {v- M) 



[v - M)+ dx 
+ dx 



< -m{n{v{t) > M)) - D2\\{v - M)+f + D2 



n(v{t)>M,z(t)>M) 



{z - M)+{v - M)+dx, 



(4.5) 



where we noticed that 



n{v(t)>M,z{t)<M) 



{z - M){v - M)+dx < 0. 
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Similarly, by taking the inner-product ( ()1.4p . {z{t) — M)_|_), where M is given in ()4.4|) . and 
through parallel steps we can get 

~\\{z - M)+f + d2 [ |V(z - M)+|2 dx < ^m{n{z(t) > M)) 

(4.6) 

- D2\\{z - M)+f + D2 [ {z- M)+iv- M)+dx, 

Jn{v(t)>M,z{t)>M) 

Sum up (|4.5p and (|4.6p and then use (|4.4p to obtain 



I {\\{v - M)+f + \\{z- M)+f ) + 2d2 iWViv - M)+f + \\Viv - M)+f ) < 



2D2 I \\{v - M)+f - 2 [ 
\ Jn{v 



{z - M)+{v - M)+ dxf + \\iz- M)+f 

(t)>M,z(t)>M) 



62 
< — e. 
- 2 



By Poincare inequality and Gronwall inequality, it follows that, for t>0,go£ Bq, 

\\(v{t)-MUf + \\{z{t)-MUf < (11(^0 - + \\{z, - M)+f ) (4.7) 

Thus there exists a time 7+(e) > Tq such that for any t > and any go £ Bq, one has 

||(r;(t) - M)+f + \\{z{t) - MUf < (4.8) 

270(2 

Symmetrically we can prove that there exists a time r_ (e) > Tq such that for any t > T_ 
and any go £ Bq, one has 

II (t;(t) + M)_ IP + II (z(t) + M)_ f < i!^, (4.9) 



Adding up (gS]) and ([O]) . we find that 

[\z{t)\-Mfdx < -"^^ 

in{\v{t)\>M) Jn{\z{t)\>M) 



[ {\v{t)\ - Mf dx + f {\z{t)\-Mfdx<^, (4.10) 

Jn(\v(t)\>M) Jn(\z(t)\>M) 7"2 



for any t > Ti = max{r+, T_} and for any go £ Bq. 

Moreover, since for any go S Bo and any T > Tq, we have 

min{\vit)\ > kM)) + m{n{\z{t)\ > kM)) < 
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there exists a sufficiently large integer k > such that for any t > Ti and go & Bq, it holds 
that 

\v{t)\'^dx+ / \z{t)\'^dx 

n{\v{t)\>kM) Jn{\z{t)\>kM) 

<2 I {\v{t)\ -Mfdx + 2M'^m{n{\v{t)\>kM)) 

Jn{\v(t)\>M) 

' (4.11) 

+ 2/ {\z{t)\ -Mfdx + 2M'^m{Q{\z{t)\>kM)) 

Jn{\z(t)\>M) 

2}?e 2M^Ko _ 2b^e 2Ko ib^e 

Therefore ()4.3p is proved with Mi = Mi(e) = kM, where M is given in ()4.4p and is the 
integer that validates (j4.1ip . and Ti = Ti(e) = max{T+,T_}. □ 

This lemma shows that the condition (jS.Sp in the item (ii) of Lemma [5] is satisfied for 
any given M > Mi(e) and any T > Ti{e), where Mi and Ti are given in Lemma [H 

Let Py : H ^ L'^{Q,)y and Pz : H ^ Lp'{Q)z be the orthogonal projections from H onto 
the T;-component space and the z-component space, respectively. The next lemma is to 
check the condition (13. 5[) and (jS.Gh for the [v, z) components in item (iii) of Lemma [H 

Lemma 7. For any given M > 0, it holds that 

K {Pv{.S{t)Bo)QQy(^t)\<M)) — > 0, as t^ oo, (4.12) 

K {Pz{S{t)Bo)^Qz{t)\<M)) — ^ 0, as t^oo, (4.13) 

in the space L'^{Q), where (»S'(t)i?o)rj(|?;(t)|<M) ^^i^d («S'(t)i?o)r2(|z(t)|<M) have been specified in 
(13771) and (IMD- 

Proof. Taking the inner-product ( (jl.2p . —Av{t)), we have 

-{vt,Av) +d2\\Av\\'^ = {u^v,Av) -b{u,Av) - D2{z - v,Av). 
By Green's formula and the homogeneous Dirichlet boundary condition, we obtain 
Id b^ d 

n^rW^'^W'^ + d2\\Avf < {u^v,Av) + Trrlluf + ^||Avf - D2{z,Av) - D2\\Vvf, 
2 dt 2a2 2 

where 

{u^v,Av) = — / n^lVvpdx — 2 / uv(\7u ■ \/v) dx 
Jn Jn 

= - \uVv + vVul"^ dx + / v^\Vu\'^dx< / v^\V 
Jn Jn Jn 



up dx. 



Consequently we get 

d r h'^ 

di 



r b'^ 

Vvf +d2\\Avf <2 / v'^\Vu\'^dx + —\\uf-2D2{\\Vvf + {z,Av)). (4.14) 
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Similarly, we can get the following inequality for the z-component, 

^\\Vzf + d2\\Azf <2 I z^\Vw\^dx + ^\\wf -2D2{\\Vzf + {v,Az)). (4.15) 
dt Jn d2 

We can also establish the inequality similar to (j4.14p but with integrals over the set 
^\v\,M = ^(1^(^)1 < ^) s-iid the inequality similar to (|4.15|) but with integrals over the set 
^\z\,M = < Then sum up the two to obtain 

'||V.||^,_ +||V.||^,_,UdJ||AH|L,,., + l|A^llE, 



b' 



|2 

\n,, 



, V™H,m + II^'II^„.m)+2 / v'\Vu\'ds + 2 I z^\Vw\''dx 
2i?2 (ll V^ll^i^i^,, + ||Vz||2 1^,^^^ + (z, A^;)n|„|,„ + {v, Az)^,,,,,,) (4.16) 



2/^2 II 2 '^2 II A II 2 2/^2 II II 2 '^2 II A ||2 



Since ^ + \\v\\^^ ^ < \\S{t)gQ\\^ < Kq due to ()4.2p . and by Poincare inequality, it 

follows that 



<^{b' + 2Dl) + 2M2 (||Vn||2 ^^i^^^ + \\Vw\\\^^^^^ , t > Tq, 50 G i^o- 
This inequality (j4.17p implies that 



(4.17) 



where 



^<rp + h, t>To, (4.18) 

= l|Vt;||^,„,,, + ||Vz||2|^l^^,^, r{t) = ^^d2, and 
/i(t) = ^{b' + 2Dl) + 2M2 (||Vn||2 ^^i^^^ + \\Vw\\l^^^^^^^ 
By (j2.17p . there exists a constant T2 = T2{Kq) > such that T2 > Tq and 

(llV^(s)ll^^^^^^^ + ||Vz(.)||2 1^,^^^) < r' {\\yv{s)f + \\Vz{s)f) ds 

2 (4.19) 
< = ^ f 1 + ^) , for t > T2, go G i?o. 

d2 V 7"2 / 
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By integrating the inequality ()2.4p on the time interval [t, t + and using ()4.19p and (|2.23p . 
we can deduce that there exists T3 = T3(Kq) > such that T3 > T2 and 



452 



<i?i + ^^^i^C. + 4Ko+(i£ + 8a2)N, t>Ts, 



(4.20) 



where Ri is the constant given in ()2.23p . 

Similarly, doing the same to (I2.10p and using ()4.19p and ()2.28p , we find that there exists 
T4 = T4(Kq) > such that T4 > T2 and 

di r' \\Vi^{s)f ds < + ^^i^^Cr + 2|1 + 2{D, - D^)!' f i^o + tt^I^^I 



d 



27^2 



< i?2 + 



2\di - d2 
di 



Cj + 2|1 + 2(Z)i - i?2)P ( Kq + ) , t>n, 



(4.21) 



where is the constant given in (j2.28p . 

From dUSD, ([OU|) and ([OT]) it follows that, for t > max{r3,r4} and any go G Bq, 



t+i 



t+i 



< 



2 

t+i 



iy{s) + ijis))-Vvis) 



ds 



||Vy(s)r ds + / ||VV'(s)r + 2 / ||V7;(s)r < Cg 



(4.22) 



where 



= 2C7 ( 1 + 



2|(ii -(i2| 



+ + i?2) + ^ (2 + |1 + 2(i^i - 1)2)1') 



+ 



— + 8a2 + 1 + 2(Di - D2f\ — 
702 / 702 



We can also assert that 
rt+l 



i+1 



^(y(s) -V'(s)) - Vz(s) 



According to (|4.22p and (|4.23p . we have 



(4.23) 



/ /i(s) ds < iM^Cs + -rib'^ + 'iDl), for t > maxfTs, r4}, go G ^o- (4.24) 

Jt "2 



Besides we have J^^^ r(s) ds < 7^2- 



Nvmk. M + .r<{C7 + ^M'Cs + -^{b' + 2Dl) e-y'- t > T,, go G Bo, 
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Finally by (j4.19p and (|4.24p and applying the uniform Gronwall inequality [34, 36] to 
(|4.18p . we obtain 

^\v\,M ' II ' ii^'|z|,M — y ' ' ^° ' ^2 

(4.25) 

where T5 = maxlTs, T4} + 1. Note that the right-hand side of ()4.25p is a uniform constant 
depending on the constant Kq in ()2.15p and the arbitrarily fixed constant M only. The 
inequality (I4.25P shows that for any given t > T^, 

Pv{S{t)Bo)fi(^\v(t)\<M) is a bounded set in Ho{^l), and 

Pz{S{t)Bo)n{\z{t)\<M) is a bounded set in H^{fl). 

Due to the compact Sobolev embedding Hq{Q,) ^ L^(il) for space dimension n < 3, it 
shows that for any given t > T^, 

Pv{S{t)Bo)fi(\v(^t)\<M) is a precompact set in L^(f]), and 

Pz{S(t)Bo)n(\z(^t)\<M) is a precompact set in L'^{0.), 

By the first property of the k- measure listed in Section 1, (j4.12p and (14.13P are proved. □ 

This lemma shows that the conditions in the item (iii) of Lemma [5] are verified for the 
t'-component in ()3.5|) and for the z-component in (13. 6p . 



5. K-Contracting Property for the {u, w) Components 

In this section, we shall check that the conditions specified in the items (ii) and (iii) of 
Lemma [5] for the {u,w) components of the coupled Brusselator equations. 

Lemma 8. For any e > 0, there exist positive constants M2 = ^2(5) and Tq = TQ{e) 
such that the u-component u{t) = u{t,x,go) and the w- component w{t) = w{t,x,go) of the 
solutions of the coupled Brusselator equations ()l.ip - ()1.4p satisfy the following estimate, 



I 

Jni 



u{t)\^dx+ \w{t)\^ dx < L2e, for t > Tq, go e Bo, (5.1) 



in{\v{t)\<M2) J n{\z{t)\<M2) 

where is a uniform positive constant. 

Proof. Recall that in our notation ^7(1^(^)1 > M) and r2(|2;(t)| > M) will be denoted by 

O^j and 17^^, respectively. Taking the L^-inner-product ( (j2.3p . v(t)) over the subset ri'^, 
we get 

~jymiH+di\\vym'+\\ymlH 

(d2 - di)y(t)A(v(t) + z(t))dx + / (v(t) + z(t) + 2a)y(t) dx 
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SO that 

|l|yWII^H+rfi||Vy(t)||^H + ||y(t)||^H 



«1 ^'m ^'m 

|2 



(5.2) 

Multiply the inequality (j5.2p by e* and then integrate it on the interval [0, t] to obtain 



|2 



+ ^\^^'*/ e^(||V.(.)||^H+||V.(«)||^H) (5.3) 
+ e-* r f4||t;(s)||2 I + 4||z(s)||2 , + Sa^lJ^Hl) ds, t > 0, 50 G Bq. 

Jo \ M M J 

On the other hand, similar to (j2.2ip . we have 



d2 / eM ||Vt;(.)||^H+||Vz(.)||^H 

, ^ " (5.4) 

1 + ^ j ft'lf^Hle* + (1 + a(t)e*) i^o, * > 0, 50 G Bq. 



< 



Let e > be arbitrarily given as in Section 4. Since (j4.4p implies < e/2 and a(t) ^ 
as t — > 00, and by (j5.4p . there exists a sufficiently large ri = Ti{e) > Tq, such that for 
t > Ti, the following inequalities hold, 

e-*||uo + V0 + W0 + zoW^i < 4Koe-^ < e, 
2\di-d2\^ /■* , , ,,,2 , „^ , ,,,2 \ , ^ 2|di -d2p A , 1 



and 

r-t 



I (l|v.(,)||^,, . ||v.(»,||^,,) < (: , _i_) 



e-' / e' ^IIKsjlP,., +4||£(s)|P|.| +8a^|!!Sl) A 



Substituting these inequalities into (j5.3p . we have 

i„l - / |y(t)|2 < Ti e, for t > n, 50 G ^0, (5-5) 

M Jn{\v{t)\>M) 
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where Fi is a uniform constant given by 

Ti l + 4a +6 (^-'^ + 27(i2 ) ( did2 7^2 

Change to Jl^j and follow the parallel steps as shown above. Then the corresponding 
inequality is also valid, 

\\y{t)C.i = [ \y{t)f dx < Ti e, for t > n, 50 e Bo. (5.6) 

^'m Jn(\z{t)\>M) 

Next taking the inner-product ( i\2.9^ .tb(t)) over the subset ri^^J, similar to (|5.2p we get 

:^llV'(t)llLi+rfillvV'(OllLi + llV'(t)ii'i., 



< ^^-7^l|V(^;(t) - I., + |1 + 2(7^1 - 2^2)ni^^(t) - ^(t)ll^i.i 

< ''^^:'^^'7 iiv.(t)ii^,i + iivz(t)ii^i, 

+ 2\l + 2iD^-D2)\^(Mt)\\li.i + \ 



(5.7) 



Multiply the inequality ()5.7p by e* and then integrate it on the interval [0,t]. Then by 
conducting similar estimates we can confirm that, for t > and go G Bq, 

||V(t)||'|.| = /, , \m\^dx<e-'\\uo + vo-wo-zof ,^^ 

+ 2|1 + 2(Z)i - D2)\'e-' f (\\v{s)\\K, + lk(s)||L.|) ds. 

Jq \ "Af ^'a/ / 

Parallel to the argument from (j5.3p through (|5.5p . we can assert that there exists a suffi- 
ciently large T2 = T2(e) > Tq, such that 



'n{|2(i)|>M) 

where r2 is a uniform constant given by 



.1 = / |^(t)|2dx < r2e, for t > T2, go £ Bo, (5.9) 



27^2 /V ^lf^2 7t^2 

Change rij^ to Q^jJ and follow the parallel steps in proving (j5.9p . Then the corresponding 
inequality is also valid: 



[ \^it)\^ dx < T2 e, for t > T2, 50 G ^o- (5.10) 

Jn{\v{t)\>M) 



A COUPLED TWO-CELL BRUSSELATOR MODEL 



23 



Finally, let M2 = max{M, Mi}, Tg = max{Ti,r2}, where M is given in ()4.4p and Mi is 
the constant in Lemma [6] and (j4.3p . Now we can combine the established (jS.Sp . (|5.10p . and 
the earlier obtained result (14.31) to conclude that 



f^(bWI>M2) 

1 



\u{t)\^dx < 2 + u 



l„l + \\v 



Ml 



{y{t) + m) 



8b' 



+ 2||f 



< \\y 



+ 



J.I +2||w 



(5.11) 



< I Ti + + — ) e, for t > Tg, 50 G ^o- 



And we can also combine (|5.6p . (j5.9p . and (j4.3p to conclude that 



/ 

Jn 



n{\z(t)\>M2) 
1 



\w{t)\''dx < 2 \\w{t) + z 



+ z 



{y{t)-m) 



+ 2|U 



< l|y(t)ll^w + 

'''■Ml 



M +2|k(t)|| 



Ml 



(5.12) 



< (Ti+T2 + —]e, for t > Tg, 50 G ^o- 



Therefore, (|5.ip is proved with 



2 ri + r2 + 



7^2 



The proof is complete. 



□ 



This lemma shows that the condition p.4p in the item (ii) of Lemma [5] is satisfied for 
any given M > M2(e) and any T > Te(e), where M2 and Tg are given in Lemma [8l 

Let Pu : H ^ L''{Q)u and Pu, '■ H L^(J7)^ be the orthogonal projections from H onto 
the n-component space and the it;-component space, respectively. The next lemma is to 
check the condition (j3.5p and (|3.6p for the {u,w) components in item (iii) of Lemma O 



Lemma 9. For any given M > 0, it holds that 

K {PuiS{t)Bo)Q(\y(t)\<M)) 
K (-Pu.(5'(t)i?o)n(|z(t)|<M)) 



0, as t ^ 00, 
• 0, as t ^ 00, 



(5.13) 
(5.14) 



in the space L^(J7), where (5'(t)i?o)n(|t)(t)|<Af) fl^t?('S'(i)-^o)r2(|z(t)|<M) have been specified in 
(13771) and (ESD. 



Proof. Taking the inner-product ( (ll.ip . — An(t))n, , j^^, we can get 
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+ (6+1) / uAudx— / aAudx — Di / {w — u)Audx 

•^^\v\,M •'^|D|,Af •^^\v\,M 

<M [ n'\Au\ dx + ^\\Au\g + U{b + l)u - a\\l 

di . D? 2 

+ ^l|A«||n|„,,M + ^ll^-^ll^.i.A^ 



f D"^ 1 
<— / n^dx + dillAull?,, , H - -ullo, , H 11(6 + l)n - all^ , . 



Note that for n < 3 the Sobolev embedding Hq{U) ^ L'^{Q) is continuous and there exists 
a uniform constant (5 > such that 



ih^n) < ^\Ml^,iny ^ ^ ^oi^)^ (5-15) 
which also holds if 0, is replaced by . Thus it follows that 

d M"^ 9D'^ 4 



< ^ {M^6^\\Vu\\^ + 4Ko{ib + l)^ + Df) +Aa^\n\) fort > Tq. 



(5.16) 



The inequality (j5.16p can be written as 



in which 



d3 

-j^ <rp + h, for t > To, (5.17) 



/3(t) = l|Vn||2^^,^^^, r(t) = ^||VH|^,,„,^^,^, and 



/i(t) = -i (i^o((fe + 1)' + Dj) + a2|f)|) . 
In view of (j4.22p . we had proved that 

rt+l 

l|Vn(s)||^l^l ds < Cs, for t>n = ma^{Ts,n} + 1, 

where Ts is the same as given in the proof of Lemma [71 Then we can apply the uniform 
Gronwall inequality to (I5.17|) to obtain 

l|V^WII%,,M ^ + ^ {Koiib + If + Dl) + a^ll^l) j exp (^^C^s) , (5.18) 
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for t>T^ and Qq ?i Bq. Inequality (|5.18p shows that 

Pu(S'(t)5o)Q(|t,(i)|<A/) is a bounded set in Hl{^l), for any t > T5, 

so that 

P„(S'(t)Po)f7(|,;(t)|<Af) is a precompact set in L'^{^), for any t > T^. 

Therefore, by the first property of the K-measure listed in Section 1, (j5.13p is proved. 

Similarly, starting from the inner-product ( (|1.3p . —l^w{t))ni^^^ and using (j4.23p . we can 
also confirm that 

/ 4 \ /2/l/f2/^2 \ 

\\vw(t)\\\^,,, - + + + + "'I^D J V^^^V ' ^^'^^^ 

for t > and 170 G i?o- Therefore, it holds that 

-fto('S'(t)-Bo)r2(|z(t)|<A/) is a precompact set in ^^(ri), for any t > T5. 
Consequently (j5.14p is proved. The proof is completed. □ 

This lemma shows that the conditions in the item (iii) of Lemma[5]have also been verified 
for the ti-component in (13. 5p and for the -w-component in (13. 6|) . 

6. The Existence of a Global Attractor and Its Finite Dimensionality 

In this section we finally prove Theorem [1] (Main Theorem) on the existence of a global 
attractor, denoted by =e/, for the coupled Brusselator semiflow {S'(t)}f>o and that £/ has 
finite Hausdorff and fractal dimensions. 

Proof of Theorem 1 . In Lemma [3l we have shown that this coupled Brusselator semiflow 
{5(t)}t>o has a bounded absorbing set Bq in H and the condition (i) in Lemma [5] is 
satisfied. 

By Lemma[6]and LemmaO we have shown that {S{t)}t>o satisfies (j3.3p and (|3.4p with 
M = max{Mi, M2} and T = max{Ti, Tg}, where Mi, Ti and M2, Tg are given in Lemma[6] 
and Lemma [8] respectively, so that the conditions (ii) in Lemma [5] is satisfied. 

By Lemma [7] and Lemma El we proved that (13. 5p and (13. 6p in the condition (iii) of 
Lemmal[5]are satisfied by {>S'(t)}t>o. Finally we apply Lemma[5]to reach the conclusion that 
there exists a global attractor ^ va. H for this coupled Brusselator semiflow {S{t)}t>o- D 

Let £^ be the global attractor of the semiflow {5(t)}t>o in H. Let qm = lini sup^^o^ Qm{t), 
where ^ 

g^(t)=sup sup [ Tr{A + F'{S{T)go))oQ^{T)dT] , (6.1) 

aoe^/ 9i£H,\\gi\\=l \t Jo J 

in which Tr i^A + F'{S{T)gQ)) is the trace of the linear operator A + F'{S{T)go), with F{g) 
being the nonlinear map in (jl.l2p . and Qm{t) stands for the orthogonal projection of the 
space H on the subspace spanned by Gi(f), • • • , Gm{t)-, with 

Gi{t) = L{S{t)go)gi, i = l,---,m. (6.2) 
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Here F'{S{T)go) is the Frechet derivative of the map F at 5'(r)goi and L{S(t)go) is the 
Frechet derivative of the map S(t) at go, with t being fixed. 

The following lemma, cf. [36, Chapter 5], will be used to show the finite upper bounds 
of the Hausdorff and fractal dimensions of this global attractor £/. 

Lemma 10. If there is an integer m such that qm < 0, then the Hausdorff dimension 
dni-i^) and the fractal dimension dpis^) of £/ satisfy 

dui'S^) < m, and dpi^/) < m max (1 + ) < 2m. (6.3) 

l<j<m-l \^ \qjn\ ) 

It can be shown that for any given t > 0, S(^) is Frechet differentiable in K and its 
Frechet derivative at g^ is the bounded linear operator L(S(t)gQ) given by 

L(5(t)5o)Go =' G{t) = {U{t), V{t), W{t), Z{t)), for any Go = (C/q, Vo, Wo, Zq) G H, 
where {U{t),V{t),W{t), Z{t)) is the strong solution of the following coupled Brusselator 



variational equation 
dU 

— = diAU + 2u{t)v{t)U + u'^{t)V -{b + 1)U + Di{W - U), (6.4) 
8V 

— = d2AV - 2u{t)v{t)U - u^{t)V + hU + D2{Z - V), (6.5) 
dW n 

— = diAW + 2w{t)z{t)W + w^{t)Z - (6 + 1)W + Di{U - W), (6.6) 
dZ 

— = d2AZ - 2w{t)z{t)W - w'^{t)Z + hW + D2{V - Z), (6.7) 

U{0) = Uo, V{0) = Vo, W{0) = Wo, Z{0) = Zo. (6.8) 

Here g{t) = {u{t),v{t),w{t), z{t)) = S{t)gQ is the solution of (jl.l2p with the initial condition 
g{0) = go- The initial value problem (16.6p ~ (|6.8p can be written as 

dC 

= (A + F'{S{t)go))G, G(0)=Go. (6.9) 



Note that the invariance of implies ^ C Bq, where Bq is the bounded absorbing set 
given in Lemma [3l Hence we have 

sup \\S{t)gof < Ko, 

Lemma 11. For the global attractor si of the coupled Brusselator semiflow {S{t)}t>Q, 
there exists a uniform constant Ki > such that 

\\Vgf<Ki, forany^Gjz/. (6.10) 

Proof. In (frT2]) . A : D{A){= U) ^ H is a positive sectorial operator and F G C^^'^^E, H). 
By (|1.13p . for any ^ ■• there is a to ^ (Oi 1/2) such that S{tQ)gQ £ E. By the solution 
theory in [34, Section 4.7], one has 

5(-)5o G C([to,cx)),S)nC7°'J((to,cx)),^)nC((to,oo),n), (6.11) 
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0,1 

were C^^^ stands for the space of Holder strongly continuous functions with exponent 1/2. 

Since S{t)^ = £/, for any g G £/ and any t > 1, there is a particular go £ £/ such that 
g = S{t)gQ. Therefore, the global attractor has the regularity that £/ C E. 

Next we prove that is a bounded set in E, so that (I6.10p holds. Suppose the contrary. 
Then there exist sequences {A''^} C (0, oo), with Ni > £, and {gi} C £/, such that 

||V<7^||>iV,, i = l,2,---. 

Let £ £/ he given such that gi = S{Ni)g^. By (|6.1ip . there is a Holder constant cq > 
such that 

0||2 ^ I Ar Co \ ^ ^ ^ ^ 1 , 1. 



||V5(t)5,ir > ^Ne--^j^, for teIe = {Ne--,N,+ 
This shows that 

||V5(r)g^°f(ir> iV^--^ — > oo, as ^ ^ oo. 
This is a contradiction to the fact that for any go G due to (I2.17|] . 

(||Vt;(s)f + ||Vz(s)f)d,<-(^i^o + ^J+^, foranyt>0, 

and, due to K22\\ and (031) . 

i-t+i 

y (||VM(s)f + ||Vu;(s)f )ds < 2C8, for any t > 0, 

where the right-hand side of the above two inequalities are uniform positive constants. 
Therefore, the conclusion holds. □ 

Theorem 2. The global attractors si for the coupled Brusselator semiflow {S'(t)}t>o has 
a finite Hausdorff dimesion and a finite fractal dimension. 

Proof. By Lemma [TOl we shall estimate Tr (^4 + F' (S{T)gQ)) o Qm{T). At any given time 
r > 0, let {(pjir) : j = 1, - ■ ■ , m} be an i/-orthonormal basis for the subspace 

Q^(r)/7 = Span{Gi(T),--- ,G,(r)}, 

where Gi{t), • • • , Gm{t) satisfy (j6.9p with the respective initial values Gi,o, • • • , Gm,o and, 
without loss of generality, assuming that Gi,o, • • • , Gm,o are linearly independent in H. By 
the Gram-Schmidt orthogonalization scheme, (pjir) = {(p^ (t) , ipj (t) , ip^ (t) , ip'j (t)) G E, for 
j = 1, • • • ,m, and ^j{T) are strongly measurable in r. Let do = mm{di,d2}. Then we 
have 

m 

TT{A + F'{S{T)go)oQ^{T) = ^ ((A(^,-(r), v,,(r)) + (F'(5(r)5o)^,(r), v^,(r))) 

(6-12) 

< -doY^ l|Vv?j(r)f + J1 + J2 + J3, 
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where 

m „ 

Ji=Y. 2u{t)v{t) {\^]{t)\'' - ^]{T)^]{r)) dx 
j=i 

+ J2 [ 2w{r)z(r) {\^'^{r)\' - v^f (T)^|(r)) dx, 

m „ 

J^=J2 (^'(^) (^](^)^|(^) - 1^1(^)1') + ^'(^) {'P'A^)^jir) - dx 

j=l •'^ 

<f2 [ {u\r)\y.]{r)\\^]ir)\+w\T)\y.^{T)\\^^{r)\)dx, 
j=i 

and 

m „ 

^3 = E / {-{b+l){\^}{r)\' + |v.|(t)P) +6(^](r)vp|(r) +(^|(r)^j(r))) dx 

m „ 

- E / (^1 {^}ir) - ^](t))' + D, (^](t) - v.|(r))') 

We can estimate each of the three terms as fohows. First, by the generalized Holder 
inequality and the Sobolev embedding Hq{Q,) L^(f2) for ra < 3, and using Lemma \TT\ 
we get 

m 

Ji <2j2Mr)h4v{r)\\L^ {y]{r)\\% + y]{r)\\L4v]{r)\\L^) 
i=i 

m 

+ 2Y,\\w{T)h44T)U. (Ilv^|(r)||i4 + ||^,^(r)||i4||v.j(r)||^4) 

(6-13) 



m m 

IIC>^,\„_I|2 ll,„V,M|2 / llv7C/,^„_ Il2||,, /,M|2 

L4 



j=i j=i 



m 



where 6 is the Sobolev embedding coefficient given in (jS.lSp . Now we apply the Garliardo- 
Nirenberg interpolation inequality, cf. [34, Theorem B.3], 

WKv < CWipCr^Mll'' , for ip G W-'^(O), (6.14) 
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provided that p,q,r > 1,0 < 6 < 1, and 

k < 9 { m I — (1 — 9) — , where n = dimil. 

p \ q J T 

Here with W^^p{Q) = L^(0), = H^{n),U'{n) = L'^{VL), and ^ = n/4 < 3/4, it 

fohows from (16.14P that 

||</',(r)||i4 <C||Vv?,(r)||t||v5,(r)||i-f =C7||V(^,(r)||i, j = !,••• ,m, (6.15) 

since = Ij where C is a uniform constant. Substitute (j6.15p into (|6.13p to obtain 

m 

Ji <45i^iC2^ ||V(/p,(r)||t. (6.16) 
i=i 

Similarly, by the generalized Holder inequality, we can get 

m m 

J2 < <5ifi^||(^,(r)||i4 < 5K^C^Y. llV'/^^-WII^- (6-17) 
i=i j=i 

Moreover, we have 

m 

i=i 

Substituting (l6lB . (leTTD and (l638|) into (|6l^ . we obtain 

m m 

Tr(A + F'(5(r)go)oQm(T) < -do ||V(/.,-(r)||2 + SJi^iC^ ^ ||V(/p,(T)||t + 6m. (6.19) 
By Young's inequality, for n < 3, we have 

m , m 

5(5i^iC^ J^||V^,(r)||t < ^^||V(^,(r)f + i^2(n)m, 
i=i i=i 
where K2{n) is a uniform constant depending only on n = dim (J7). Hence, 

, m 

TV (A + F'(5(t)5o) o Qm{T) < -y E + (^2(n) + 6) m, r > 0, 50 G =2^- 

i=i 

According to the generalized Sobolev-Lieb-Thirring inequality [36, Appendix, Corollary 
4.1], since {(^i(r), • • • , V'm(''")} is an orthonormal set in H, so there exists a uniform constant 
> only depending on the shape and dimension of O, such that 

m II 2 

m 1 



^||Vy,,(r)|p>i^3^. (6.20) 

Therefore, we end up with 

Ti {A + F' {S{T)gQ) o Q^{t) < -^^m^+l + {K2{n)+b)m, T>0,go£s^. (6.21) 

2\n\n 
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Then we can conclude that 

ft 



q^{t)=SMV sup / T¥(^ + F'(5(r)5o)) oQ^(r)dr) 

9oe.e/ gr£H,\\gi\\=l JQ J 



< — — — f-m"'^'''" + (K2(n) + b) m, for any t > 0, 
2101 



(6.22) 



so that 



hmsup Qmit) < - -^^m^+n + (i^2(n) + 6) m < 0, (6.23) 



if the integer m satisfies the following condition, 



,'2(K2(n) + 6)\"/^^, , , 

m - 1 < ^ \' M |r?|<m. (6.24) 



4^3 , 

According to Lemma \W\ we have shown that the Hausdorff dimension and the fractal 
dimension of the global attractor £/ are finite and their upper bounds are given by 

dni-s^) < m and dp{.sif) < 2m, 

where ra satisfies (I6.24p . □ 
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